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We present a method to perform fault-tolerant single-qubit gate operations using Landau-Zener 
tunneling. In a single Landau-Zener pulse, the qubit transition frequency is varied in time so that 
it passes through the frequency of the radiation field. We show that a simple three-pulse sequence 
allows eliminating errors in the gate up to the third order in errors in the qubit energies or the 
radiation frequency. 
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I. INTRODUCTION 

In many proposed implementations of a quantum com- 
puter (QC) single-qubit operations are performed by ap- 
plying pulses of radiation. The pulses cause resonant 
transitions between qubit states, that is between the 
states of the system that comprises a qubit. The opera- 
tion is determined by the pulse amplitude and duration. 
In many proposals, particularly in the proposed scalable 
condensed-matter based systems jlj, control pulses will 
be applied globally, to many qubits at a time. A target 
qubit is chosen by tuning it in resonance with the radi- 
ation. The corresponding gate operations invariably in- 
volve errors which come from the underlying errors in the 
frequency, amplitude, and length of the radiation pulse 
as well as in the qubit tuning. 

Improving the accuracy of quantum gates and reducing 
their sensitivity to errors from different sources is critical 
for a successful operation of a QC. Much progress has 
been made recently in this direction by using radiation 
pulses of special shape and composite radiation pulses 
[21 • In the analysis or resonant pulse shape it is usu- 
ally assumed that the qubit transition frequency is held 
constant during the pulse. 

An alternative approach to single-qubit operations is 
based on Landau-Zener tunneling (LZT) [^Q. In this 
approach the qubit transition frequency u>o{t) is swept 
through the frequency of the resonant field cop 0- The 
change of the qubit state depends on the field strength 
and the speed at which uJo{t) is changed when it goes 
through resonance The LZT can be used also for a 
two-qubit operation in which qubit frequencies are swept 
past each other leading to excitation swap 0, S ■ 

In the present paper we study the robustness of the 
LZT-based gate operations. We develop a simple pulse 
sequence that is extremely stable against errors in the 
qubit transition frequency or equivalently, the radiation 
frequency. Such errors come from various sources. An 
example is provided by systems where the qubit-qubit 
interaction is not turned off, and therefore the transition 



energy of a qubit depends on the state of other qubits. 
Much effort has been put into developing rneans for cor- 
recting them using active control 0, llfl UM ■ 

An advantageous feature of LZT is that the change 
of the qubit state populations depends on the radia- 
tion amplitude and the speed of the transition frequency 
change luq, but not on the exact instant of time when 
the frequency coincides with the radiation frequency, 
u!o{t) = iup. However, the change of the phase differ- 
ence between the states depends on this time. Therefore 
an error in luq or up leads to an error in the phase dif- 
ference, i.e., a phase error. This error has two parts: 
one comes from the phase accumulation before crossing 
the resonant frequency, and the other after the crossing. 
Clearly, they have opposite signs. 

A natural way of reducing a phase error is to make 
the system accumulate the appropriate opposite in sign 
phases before and after the "working" pulse. To do this, 
we first apply a strong radiation pulse that swaps the 
states, which can be done with exponentially high efR- 
ciency using LZT. Then we apply the "working" pulse, 
and then another swapping pulse. The swapping pulses 
effectively change the sign of the accumulated phase. As 
we show, by adjusting their parameters we can compen- 
sate phase errors with a high precision. 

In Sec. II below we give the scattering matrix for LZT 
in a modified adiabatic basis which turns out to be ad- 
vantageous compared to the computational basis. The 
scattering matrix describes the quantum gate. In Sec. Ill 
it is presented in more conventional for quantum compu- 
tation terms of the qubit rotation matrix. In Sec. IV, 
which is the central part of the paper, we propose a sim- 
ple composite Landau-Zener (LZ) pulse and demonstrate 
that it efficiently compensates energy offset errors even 
where these errors are not small. Sec. V contains con- 
cluding remarks. 



II. LANDAU-ZENER TRANSFORMATION IN 
THE MODIFIED ADIABATIC BASIS 
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A simple implementation of the LZ gate is as follows. 
The amplitude of the radiation pulse is held fixed, while 
the difference between the qubit transition frequency and 
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the radiation frequency 



A{t) — Lup — u>o(t) 



(1) 



is swept through zero. If ujQ{t) is varied slowly compared 
to ojF, i.e., \uJo\ -C ujp, the qubit dynamics can be de- 
scribed in the rotating wave approximation, with Hamil- 
tonian 



(2) 



Here, 7 is the matrix element of the radiation-induced 
interstate transition. The Hamiltonian H is written in 
the so-called computational basis, with wave functions 

|0)= (^J) and|l)= (^J 

We assume that well before and after the frequency 
crossing the values of |A| largely exceed 7 and A slowly 
varies in time, |A/A^| ^ 1. Then the wave functions of 
the system are well described by the adiabatic approxi- 
mation, i.e., by the instantaneous eigenfunctions of the 
Hamiltonian ((^J, 



1^0) 



cos(6'/2) 
sin(6'/2) 



- sin(6'/2) 
cos(6l/2) 
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where A = A(<) and (— l)"ii^sgnA is the adiabatic en- 
ergy of the states = The adiabatic approxi- 
mation for E and 9 is accurate to 7^A/A'^ and 7A/A'^, 
respectively. 

In contrast to the standard adiabatic approximation, 
we chose the states |V'o,i) and their energies in such a way 
that IV'o) and IV'i) go over into |0} and |1), respectively, 
for |A|/7 — !■ cx). As a result 6 is discontinuous as a func- 
tion of A for A = 0, but the adiabatic approximation 
does not apply for such A anyway. 

For the future analysis it is convenient to introduce the 
Pauli matrices X, Y, Z in the basis Q, with 

Z\4,n) = {l-2n)\^/M), X\^n) = \^i^n) (ri-0,1), 

and Y — iXZ. In these notations, the operator of the 
adiabatic time evolution U{tf, ti) = Tcxp[— i j^' dtH(t)] 
has the form 



U(tf, ti) = exp 



Eit) dt 



(4) 



with sgnA = sgnA(ti) = sgnA(ij) [the sign of A(i) is 
not changed in the range where Eq. applies] . 

The LZ transition can be thought of as occurringbe- 
tween the states Following the standard scheme 
we take two values Ai 2 of A(i) such that they have op- 
posite signs, A1A2 < 0. We choose |Ai 2I sufficiently 
large, so that the adiabatic approximation Q applies for 



A(ti) = Ai,i = 1,2. At the same time, |Ai^2| are suffi- 
ciently small, so that A(t) can be assumed to be a linear 
function of time between Ai and A2, 



A(0 « -77(i-ic), 7] = -l\{t, 



(5) 



where the crossing time tc is given by the condition 
A(ic) — 0. The adiabaticity for t = ^1.2 requires that 
Ai^2| ^ 7, ?7^^^. We will consider the LZ transition first 
for the case Ai > and A2 < 0, when 77 > 0. 

The modified adiabatic basis J^Jl is advantageous, be- 
cause in this basis the transition matrix S has a par- 
ticularly simple form. For A(<) of the form ^ the er- 
ror in 5* is determined by the accuracy of the adiabatic 
approximation itself and is of order ^l'rf'\t\_-i ~ tcP, in 
contrast to the computational basis, where the error is 
^ 0(7/771^1.2 — tc\). This latter error is comparatively 
large for the values of 7/IA1 2I of interest for quantum 
computing. It leads to the well-known oscillations of the 
transition amplitude with increasing |A| 6], whereas in 
the basis © such oscillations do not arise, see Fig^ 



1.0 



TO 

o 



q 
'.4—' 
w 

03 



0.5 



0.0 



-10 




FIG. 1: (Color online) . Landau- Zener transitions |i/>o} li^i) 
and |0) |1) in the modified adiabatic basis and in 
the computational basis for linear A{t) Solid and dot- 

ted lines show the squared amplitude of the initially empty 
states \tpi) and respectively. The data for the state \tpi) 
close to A{t) = are not shown, since the adiabatic ap- 
proximation does not apply for small |A|. The lines refer 
to g — 1, 0.47, 0.33, and 0.21, in the order of decreasing tran- 
sition probability for {t — tc)ri^^^ = 10. As long as A{t) is 
large and negative, the system stays in the initially occupied 
adiabatic state \tpo), and therefore the solid curves for dif- 
ferent g cannot be resolved for (t — tc)r]^^^ < —1. For large 
{t~tc)rj^^'^ the solid lines quickly approach the Landau-Zener 
probabilities 1 — e 



The energy detuning |Ai_2| cannot be made too large, 
because this would make the gate operation long. If we 
characterize the overall error of the adiabatic approxima- 
tion as the sum X)i=i 2 ^/v'^l^i ~^cP and impose the con- 
dition that the overall duration of the operation t2 — ii 
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be minimal, we see that the error is minimized when 
the pulses A{t) are symmetrical, <2 — = tc — ti, i-e., 
|Ai| = |A2|. 



The matrix S{t2,ti) = S' in the basis (PJ can be ob- 
tained using the parabolic cylinder functions that solve 
the Schrodinger equation with the Hamiltonian (jSJ, lO, 



exp [-7rg2 + i[ip^ - i^^)] 

i{ipi + 1P2)] 



gr(-ig2) ^-^y [ 2 .9 + * 4 



(27r)^/" 

exp [- 



(6) 



where T{x) is the gamma function. 

The dimensionless coupling parameter g = 7/|?y|^^^ in 
Eq. ^ is the major parameter of the theory, it deter- 
mines the amplitude of the {ipn) — IV'i-n) transition. 
The phases (pi^2 are 



4|ry| 



'In 



1/2 



2A? 



1,2. (7) 



Here we have disregarded the higher order terms in 
|Ai.2|~^. The constants in (pi^2 are chosen so as to 
match the corresponding constants in the parabolic cylin- 
der functions |ia |. 

The matrix S for a transition from the initial state 
with Ai < to the final state with A2 > is given by 
the transposed matrix jnj in which the phases ipi and 
ip2 are interchanged. In this case 77 < in Eq. lO; the 
expressions for <pi_2 and g do not change. 



The rotation angles a are given by the expressions 



<P = -2ipo + argViig^) + 
a — 2cos~^ [exp(— 7rg^)] 



(11) 



A minor modification of these equations allows using 
them also for the case Ai < < A2 when the frequency 
difference is increased in time in order to bring the states 
in resonance. It was explained below Eq. |(7J how to 
relate the matrix S in this case to the matrix S for Ai > 
> A2. Following this prescription we obtain 

S{t2,ti) = [20(i2) - *] RAa)Rz h2(/.(ti) + $] . 

(12) 

In the rotation matrix representation, the only difference 
from the S matrix from the case of decreasing A(t) is that 
$ and (j){ti,2) change signs. Eqs. lp^- (fT^ express the LZ 



III. ROTATION MATRIX REPRESENTATION 

The LZ transition can be conveniently described us- 
ing the standard language of gate operations in quantum 
computing. To do this we express the transition ma- 
trix in terms of the operators Rx{Q) = exp(— i6'X/2) and 
Rz{6) — exp{~i6Z/2) of rotation about x and z axes in 
the basis © . The rotation matrices can be written using 
the "adiabatic" phases 4>{ti) that accumulate between the 
time ti and the time tc at which the levels would cross in 
the absence of coupling. From Eq. {Tj) 

ipi ^ (f)iti) + ipa (i = l,2), ti<tc<t2, 

1 



Edt 



V(lng^-l), (8) 



where we have disregarded corrections oc |Ai^2| ^, in 
agreement with the approximations made in obtaining 
Eq. ©■ 

For the case Ai > > A2 the dependence of the tran- 
sition matrix S" ((HJ on the phases (?!'(ti,2) has the form 



5(i2, ii) = Rz S'R, [2(/)(ti)] . (9) 

A direct calculation shows that the matrix S" is 

5' = i?,($)i?,(a)i?,(-$). (10) 




FIG. 2; (Color online). The rotation angles a (solid line) 
and $ (dashed line) in the rotation-matrix representation of 
the Landau-Zener gate operation as functions of the con- 
trol parameter g. The 7r/2 gate, a = n/2, requires g = 
[In2/(27r)]-^''2 ^ g.SS. 

transition matrix in the form of rotation operators in the 
basis of the modified adiabatic states |V'o) and lipi) (|3Jl. 
For strong coupling, exp(— Trg^) <C 1, the rotation angle 
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a approaches tt, which corresponds to a population swap 
between the adiabatic states. It is well known from the 
LZ theory 0,13 that the swap operation is exponentially 
efficient, tt — a « 2exp(— Trg^) for large g. In the opposite 
limit of weak coupling, g ^ 1, the change of the state 
populations is small, a « {SnY^^g. In addition to the 
change of state populations there is also a phase shift 
that accumulates during an operation. The dependence 
of the angles a and (f> on the coupling parameter g is 
shown in Fig [3 



IV. COMPOSITE LANDAU-ZENER PULSES 

For many models of quantum computers an important 
source of errors are errors in qubit transition frequencies 
luq- They may be induced by a low- frequency external 
noise that modulates the interlevel distance. They may 
also emerge from errors in the control of the qubit-qubit 
interaction: if the interaction is not fully turned off be- 
tween operations, the interlevel distance is a function of 
the state of other qubits. In addition there are systems 
where the interaction is not turned off at all, like in liquid 
state NMR-based QC's. In all these systems it is impor- 
tant to be able to perform single-qubit gate operations 
that would be insensitive to the state of other qubits. 

The rotation-operator representation suggests a way to 
develop fault tolerant composite LZ pulses with respect 
to errors in the qubit transition frequency loq and in the 
radiation frequency ojp- We will assume that there is 
a constant error e in the frequency difference A(t) = 
ujF—<jJo{t)-i but that no other errors occur during the gate 
operation. From Eq. ©, the renormalization A(t) — > 
A{t)+e translates into the change of the adiabatic energy 
E and the crossing time tc, with tc ^ tc+e/rj. As a result 
the phases 4>{ti^2) as given by Eq. ^ are incremented by 



to second order in e. 



|Afe)| 2 
8\v\E{U) ' 



z = l,2, (13) 



A. Error compensation with 7r-pulses 

A simple and robust method of compensating errors 
in (p is based on a composite pulse that consists of the 
desired pulse sandwiched between two auxiliary pulses. 
Using TT-pulses in which A{t) is linear in t, as shown 
in Fig. it is possible to eliminate errors of first and 
second order in e. The goal is to compensate the factors 
Rz[±26(j){ti^2)] in the ^'-matrix We note that all 

other factors in S are not changed by the energy change e, 
which is one of the major advantageous features of the LZ 
gate operation. A 7r-pulse is obtained if exp{Trg^) ^ 1, 
which is met already for not too large g: for example, 
exp(-7rg2) < 10-5 for g > 1.92. 



Disregarding corrections ~ exp(— tt^^) we can write the 
^-matrix for the 7r-pulse as 

S^it',t) « -zXi?,[20^(i) + 20,(i')-2$] 

= -iR,[-2Mt) " ^Mt') + 2$]X, (14) 

where t, t' are the initial and final times, and the sub- 
script TT indicates that the corresponding quantities refer 
to a TT-pulse. We assume that A(t) > > A(i'). 

The overall gate operation is now performed by a com- 
posite pulse 



Scit2, t'l) — S-!r{t2,t2)S{t2; ti)S-!r{ti, t[] 



(15) 



In writing this expression we assumed that the system 
is switched instantaneously between the states that cor- 
respond to the end (beginning) of the correcting pulse 
and the beginning (end) of the working pulse S{t2,ti). 
The overall composite pulse is shown in Fig. 13 The first 




FIG. 3; (Color online). An idealized composite pulse. The 
first and third pulses are 7r-pulses, the pulse in the middle 
performs the required gate operation. The overall pulse com- 
pensates errors in the qubit energy to 3rd order. 

and the second 7r-pulses correct the errors Set) ^ in the 
phases (j){ti) and (j){t2), respectively. We show how it 
works for (j>{t2)- From Eqs. (|14|l . the error in (j){t2) 
will be compensated if 



5^^{Q + 6^^{t2) - 5^{t2) = 0. 

To second order in e, the errors Scj) here are given by 
Eq. (|13|l with appropriate U. The total error will be 
equal to zero provided 



E^it'^) E{t2) E^{t2 



|A.(ty| \A{t2) 



Ax(t2) 



7y7ri?7r(4) VEit2) 1l7rErrit2 



0, 



= 0. 



(16) 



Equations (|16|l are simplified if we keep only the lowest 
order terms with respect to in which case E(ti) « 
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|A(ti)|/2 both for the working and the correcting pulse. 
This gives 

= 2r/, |A,(4)| -2|A(t2)| -A,(t2) =0. (17) 

An immediate consequence of Eq. H17|l is that the cou- 
phng constant ^.^ for the 7r-pulse should exceed the 
value of 7 for the working pulse, because g-n > g and 
r]T^ > rj. Another consequence is that the 7r-pulse am- 
plitude should exceed that of the working pulse. If we 
choose A^ so that the error of the adiabatic approxima- 
tion in the 7r-pulse does not exceed that of the work- 
ing pulse, 7^?7^/|A^p < 777/|A(t2)P, we obtain from 
Eq. lO the condition A^(t2) > \A{t2)\2^^^{g^/gY^^ . 

We note that the correcting pulse is asymmetric, with 
|A„(i^)| > A„(i2),2|A(t2)|, as shown in Fig. 01 Another 
important comment is that the proposed simple single 
pulse docs not allow us to correct errors of higher or- 
der in e. It is straightforward to see that the equation 
for A^(t2), A7r(t2) that follows from the condition that 
the error ~ vanishes is incompatible with Eqs. (|17|l . 
However, the terms cx (e/7)'^ contain a small factor 
f;^7/A(ti 2)'^ <SC 1. The higher-order terms in 2/7 con- 
tain higher powers of the parameter 7/A(ti.2). This is 
why compensating errors only up to the second order in 
£ turns out efficient. 

The analysis of the first correcting 7r-pulse, S'7r(ti, t'^), 
is similar to that given above. The amplitude of this pulse 
also exceeds the amplitude of the working pulse. The 
duration of the correcting pulses is close to the duration 
of the working pulse, for (7 ~ 1 and gj^ > 2. 

The pulse sequence (|15|) is written assuming that the 
radiation is switched off between the pulses and that the 
switching between the working and correcting pulses is 
instantaneous. A generalization to a more realistic case 
of a nonzero switching time is straightforward. The time 
evolution between the pulses can be described by extra 
terms in the phases (f>Triti),4'Tr{t2), leading to the appro- 
priate modification of the equations for error compensa- 
tion Hlt)|) . The analysis can be also extended to the case 
where A(t) is a nonlinear function of time and the cou- 
pling g depends on time. This extension requires numeri- 
cal analysis; we have found for several types of A{t), g{t) 
that good error correction can still be achieved with a 
three-pulse sequence. 

B. Mciximal error of the three-pulse sequence 

In order to demonstrate the error correction we will 
consider single working pulses S{t2,ti) with the overall 
phases $ — 20(ti^2) =0(mod27r) in the absence of er- 
rors, which we will denote by S^^\t2,ti). Such pulses 
describe transformations between the states Q with no 
extra phase, that is pure X rotations. We will also 
choose the correcting 7r-pulses with the overall phase 
2(j)^{t) + 2(j){t') - 2$ = (mod27r) in the absence of er- 
rors, with i, t' being t^, ti and t2, t'2 for the first and sec- 
ond pulse, respectively. Then in the absence of errors the 



overall gate is either not affected by the correcting pulse 
or its sign is changed. 

The restriction on the phases provides extra con- 
straints on the parameters of the correcting pulses. 
First of all, it " discretizes" the total duration of the 
pulses. For the correcting pulses we still have a choice 
of A^(t2) and A7r(ii). They will be chosen maximally 
close to \A{t2)\2^/^gJgy/^ and A{h)2^/^gJ gy/^ re- 
spectively, in order to minimize the error of the adiabatic 
approximation JJJl and to minimize the overall pulse du- 
ration. 

We will characterize the gate error £ by the spectral 
norm of the difference of the operator S in the presence 
of errors and the "ideal" gate operator S^"\ 

£ = \\S - S^'^\2. (18) 

Here, ||^||2 is the square root of the maximal eigen- 
value of the operator A'^A. For uncorrected pulses S = 
S{t2,ti), whereas for corrected pulses S = Sc{t'2,t[). For 
simple symmetric composite pulses described below, the 
overall sign of the composite pulse is opposite to that of 
the original pulse in the absence of errors. In this case 
we set S = -Scit'2,t[) in Eq. l(TH|) . 
For uncorrected pulses we have 

£ = 2^/^|l - na;ina:2 - nyiny2C0sa\^/'^, (19) 

where = {cos[6(t){ti)],sm[S(f){ti)]) is an auxiliary 2D 
unit vector {i — 1,2). Eq. H18|) applies also in the case 
of corrected pulses, but now we have to replace in the 
definition of the ni vector 

6^ih) 6^ih) - (50,(ti) - Scj^M- (20) 

A similar replacement must be done in the definition of 
the vector n2. 

For small phase errors |<5'/>(ii,2)| the function £ for un- 
corrected pulses is linear in the error. In particular, to 
first order in e for a symmetric pulse we have |(50(ti)| ~ 
\S(j){t2)\ ~ eE{ti)/ri, and £ « 2\e\E{ti)Tj-^ sm{a/2). In 
contrast, by applying the same arguments to a corrected 
pulse, we see that the gate error is tx e'^. As noted above, 
the terms oc e"^ and higher-order terms in e contain a 
small factor. They become very small already for not 
too small e. 

To illustrate how the composite pulse works we com- 
pare in Fig. 2] the error of an uncorrected LZ gate with 
the gate error of the composite pulse. The data re- 
fer to different values of g of the working pulse; the 
corresponding values of a are given in Fig. |21 We 
used A{ti) = — A(t2) « lOrj^/^ [the precise value of 
A(ii,2) was adjusted to make S{t2,ti) an X-rotation, 
S(t2,ti) = Rx{a)]. The compensating 7r-pulses where 
modeled by pulses with = 3. Based on the arguments 
provided at the end of Sec. IV A, we took A^(i2) ~ 
|A(t2)|2i/2(g^/g)i/3^ whereas A^(t^) was found from 
Eq. Iinil; we used A^{t[) = -A^(i^) and A^(ti) = 
-A^{t2). 
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It is seen from Fig. 0] that the proposed composite 
pulses are extremely efficient for compensating gate er- 
rors. Even for the energy error e = 7, where the error of 
an uncorrected pulse is close to 1, for the composite pulse 
E < 10~'^. For g < 1 the error of the single pulse scales as 
£, whereas the error of the composite pulse scales as e'^, 
in agreement with the theory. For large 5, when the gate 
is almost an X-gate (7r-pulse), in the case of symmetric 
pulses that we discuss, the coefficients at the terms oc e 
and oc become small; they become equal to zero for 
a — Tx. Therefore for large g and for not too small e the 
errors of single and composite pulses scale as and ^ 
respectively. On the other hand, for 2/7 close to 1, errors 
of the composite pulses with larger g are larger than for 
smaller g. This is because the calculations in Fig. 0] refer 
to the same I^jr^l"^ ^ in which case the errors oc e'^, e'^ are 
proportional to g^ . 




sly 

FIG. 4: (Color online) . Gate errors £ for Landau- Zener pulses 
as a function of the frequency error e. The upper and lower 
curves refer to the single LZ pulse and the composite pulse, 
respectively. The dash-dotted, dotted, solid, and dashed lines 
show £ for g = 2, 1.2, 1, and 0.3. 



V. CONCLUSIONS 

In this paper we have developed a theory of quantum 
gates based on LZ pulses. In these pulses the control 
dc field is varied in such a way that the qubit frequency 
passes through the frequency of the external radiation 
field. In the adiabatic basis an LZ gate can be expressed 
in a simple explicit form in terms of rotation matrices. 
Our central result is that already a sequence of three 
LZ pulses can be made fault-tolerant. The error of the 
corresponding composite pulse £ scales with the error e in 
the qubit energy or radiation frequency at least as e^. In 
addition, the coefficient at has an extra parametrically 



small factor. The duration of the 3-pulse sequence is 
about 4 times the duration of the single pulse, for the 
parameters that we used. 

Fault tolerance of LZ gates is partly due to the change 
of state populations being independent of precise fre- 
quency tuning. In particular, LZ tunneling makes it pos- 
sible to implement simple 7r-pulses with an exponentially 
small error in the state population. 

The approach developed here can be easily generalized 
to more realistic smooth pulses, as mentioned above. It 
can be applied also to two-qubit gate operations in which 
the frequencies of interacting qubits are swept past each 
other, leading to excitation transfer Such operations 
are complementary to two-qubit phase gates and require 
a different qubit-qubit interaction. 

LZ pulses provide an alternative to control pulses 
where qubits stay in resonance with radiation for a spec- 
ified time 0- In this more conventional approach it is 
often presumed that qubits are addressed individually by 
tuning their frequencies. In contrast to this technique, LZ 
pulses do not require stabilizing the frequency at a fixed 
value during the operation. As a consequence, calibration 
of LZ pulses is also different, which may be advantageous 
for some app lications, in particular in charge-based sys- 
tems H Ell. The explicit expressions discussed above 
require that the qubit transition frequency vary linearly 
with time, but the linearity is needed only for a short 
time when the qubit and radiation frequencies are close 
to each other, as seen from FigQ] which should not be 
too difficult to achieve. 

For pulses based on resonant tuning for a fixed time, 
much effort has been put into developing fault-tolerant 
pulse sequences, see Ref. 14 and papers cited therein. In 
particular, for energy offset errors it has been shown that 
a three-pulse sequence can reduce the error to S ^ ^3 
(the fidelity F evaluated in Ref. ^15, is related to £ dis- 
cussed in Ref. and in this paper by the expression 
1 — F oc f ^ for small £; therefore an error £ ^ cor- 
responds to the estimate 1 ~ F ^ e^). This error is 
parametrically larger, for small e, than the error of the 
three-pulse sequence proposed here, £ oc . We note 
that, with two correcting pulses of a more complicated 
form, it is possible to eliminate errors of higher order in 
e. 

It follows from the results of this paper that fault- 
tolerant LZ gates can be implemented using the standard 
repertoire of control techniques and may provide a viable 
alternative to the conventional single qubit gates. 



Acknowledgments 

This work was supported in part by the NSF through 
grant ITR-0085922 and by the Institute for Quantum 
Sciences at Michigan State University. 



7 



[1] B. E. Kane, Nature 393, 133 (1998). 

[2] L. M. K. Vandersypen and I. L. Chuang, Rev. Mod. Phys. 

76, 1037 (2004). 
[3] L. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 

[4] C. Zencr, Proc. R. Soc. London, Scr. A 137, 696 (1932). 
[5] M. Dykmaii and P. Platzman, Quantum Inf. Coniput. 1, 
102 (2001). 

[6] V. Benza and G. Strini, Fortschr. Phys. 51, 14 (2003). 
[7] P. Platzman and M. L Dykman, Science 284, 1967 
(1999). 

[8] K. Saito and Y. Kayanuma, Phys. Rev. B 70, 201304 

(2004). 

[9] L. Viola, J. Mod. Opt. 51, 2357 (2004). 
[10] P. Facchi, S. Tasaki, S. Pascazio, H. Nakazato, A. Tokuse, 



and D. A. Lidar, Phys. Rev. A 71, 022302 (2005). 
[11] P. Sengupta and L. P. Pryadko, Phys. Rev. Lett. 95, 
037202 (2005). 

[12] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, 
Series, and Products, 6th edition (San Diego, CA: Aca- 
demic Press, 2000). 

[13] V. N. Smclyanskiy, A. G. Pctukhov, and V. V. Osipov, 
Phys. Rev. B 72, 081304 (2005). 

[14] K. R. Brown, A. W. Harrow, and I. L. Chuang, Phys. 
Rev. A 70, 052318 (2004). 

[15] H. K. Cummins, G. Llewellyn, and J. A. Jones, Phys. 
Rev. A 67, 042308 (2003). 



